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Abstract. Two-proton and two-deuteron correlations have been studied in the target fragmentation region of cen- 
tral Pb+Pb collisions at 158 A GeV. Protons and deuterons were measured with the Plastic Ball spectrometer of the 
WA98 experiment at the CERN SPS. The results of one-dimensional and multi-dimensional analyses using both the 
Bertsch-Pratt and Yano-Koonin-Podgoretsky parameterizations of the two-particle correlation functions are presented. 
The proton source exhibits a volume emission, while the deuteron source, with small outward radius, appears opaque. 
Both proton and deuteron sources have cross-terms i? 2 ,, and longitudinal velocities (3 consistent with zero, indicat- 
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ing a boost-invariant expansion. The invariant radius parameter R follows an approximate A / \fm scaling while the 
longitudinal and transverse radii, Rl and Rt, scale approximately as A/y/mr with A « 3 fm GeV 1/2 in both cases. 
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1 Introduction 

In this paper we report on an analysis of quantum-statistical and 
final state induced correlations [ 1 1 of proton and deuteron pairs 
emitted from 158 A GeV central Pb+Pb collisions. The study of 
these two particle species in the same system is of special inter- 
est as one would expect that deuterons, due to their loose struc- 
ture and small binding energy, will survive only in the late, low 
density environment when scatterings are rare. The generally 
accepted mechanism of deuteron production is by coalescence 
of a proton and neutron H2I3I41 at the time of freeze-out. The 
deuteron measurement thus provides an unique means to mea- 
sure the geometry of the source in the late time interval of the 
expansion. Experimental results on Fermi-Dirac correlations 
for protons are not as abundant as those on Bose-Einstein cor- 
relations for mesons, and are very scarce for deuteron correla- 
tions. A remarkable feature of the existing proton and deuteron 
correlation data in hadron-nucleus and nucleus-nucleus colli- 
sions, spanning over four orders of magnitude of incident en- 
ergy from 0.02 to 450 A GeV is the small variation of the ex- 



and made use of the Midrapidity Calorimeter(MIRAC) |23|, 
the Zero Degree Calorimeter(ZDC) 1241 . and the Plastic Ball 
spectrometer ll25ll . 

The ZDC registers energy emitted along the beam direction 
in the 3° forward cone. The MIRAC measures the total trans- 
verse energy in the pseudorapidity region of 3.5 < r\ < 5.5. 
It is a sampling calorimeter that consists of a lead-scintillator 
electromagnetic section, followed by an iron-scintillator hadronic 
section. MIRAC plays the central role in the WA98 minimum 
bias trigger where the measured transverse energy Ex is re- 
quired to be above a minimum threshold. MIRAC is used to 
classify the centrality of each event. The analysis presented 
here has been performed on the 10% most central events of the 
measured minimum bias cross section of sa 6400 mb. About 3 
millions central Pb+Pb events were analyzed with an average 
of 6.1 identified protons and 2.9 identified deuterons per event. 
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2 Experimental setup 
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Fig. 1. Experimental setup of the WA98 experiment for the 1996 run. 
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The CERN SPS experiment WA98 (for details see G2 and 
references therein) is a general purpose apparatus that consists 
of large acceptance photon and hadron spectrometers, detectors 
for charged particle and photon multiplicity measurements, and 
calorimeters for transverse and forward energy measurements. 
The layout of the WA98 experimental setup for the 1996 SPS 
run period is shown in Fig.[T] The results presented in this re- 
port were obtained from an analysis of data taken in 1996 with 
the 158 A GeV 208 Pb beam on a 239 mg/cm 2 208 Pb target 



The identification and momentum measurement of protons 
and deuterons has been made using the Plastic Ball spectrome- 
ter 11251 . The Plastic Ball consists of 655 detector modules sur- 
rounding the target over the range of polar angles d from 30° 
to 160° with full azimuthal acceptance. Each of the modules 
consists of a AE and E section. The AE section is a CaF2(Eu) 
crystal scintillator and the E section is a plastic scintillator. The 
light emission of the plastic scintillator is approximately 100 
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times faster than that of the CaF2(Eu), so signals from both 
scintillators can be read out by a single photomultiplier with 
subsequent separation of the signals by pulse shape analysis. 
The AE—E identification technique was used and is capable of 
identifying pions and charged fragments up to the helium iso- 
topes. Because of the relatively low yield of pions as compared 
to baryons, an additional positive pion identification is made 



by means of its sequential decay 7r + — > 



The thickness of the AE crystal scintillator was chosen to be 
4 mm with 35.6 cm as the length of the E plastic scintillator. 
This assures that protons up to 200 MeV are fully stopped in 
the S-counter and provide a complete AE — E signal. Addi- 
tionally, the mis-identification of punch through deuterons do 
not disturb the proton spectra below 200 MeV. 

The WA98 experimental setup has been implemented in the 
GEANT Monte Carlo simulation package [26|. This has been 
used to study the response of the Plastic Ball in order to extract 
the particle identification efficiency and the kinetic energy reso- 
lution for different particle species, as well as two-particle vari- 
ables and takes into account the granular structure of the Plastic 
Ball. GEANT transports the particles through the experimental 
setup taking into account the geometrical material boundaries 
and particle interactions. The particles are traced through the 
experimental setup where they may interact with the material 
they encounter on their way to the detector. When a particle 
reaches the detector the output signal of the Plastic Ball mod- 
ule is simulated. The particle distributions used as input to the 
GEANT simulations were adjusted to reproduce the measured 
momentum and angular distributions of particles. The Plastic 
Ball acceptance is illustrated in Fig. [2] as the distribution of 
protons and deuterons that pass the acceptance cuts. 



3 One-dimensional analysis 

The one-dimensional experimental correlation function was con- 
structed as: 

Yl2(q) 



C(q) = N 



(1) 



where q = |(qi— q2) is the two-particle three-momentum dif- 
ference in the laboratory frame, N is a normalization constant, 
and the numerator is the coincidence yield while the denomi- 
nator is the uncorrected background yield. In this analysis, the 
event-mixing technique was used to construct the background 
yields, using the same central event data set as used to con- 
struct the numerator. To insure that only particles from compa- 
rable events were mixed, the central event sample was further 
subdivided into 8 subsets according to the measured transverse 
energy. Particles were picked randomly from different events 
belonging to the same centrality selection with the additional 
condition to have the same multiplicities of the studied particle 
and all particles registered in the Plastic Ball. The number of 
mixed background pairs was chosen to be 10 times larger than 
the number of coincident pairs to insure a small statistical error 
contribution from the background pair yield, Y^ 2 . 

The measurement resolution of the variable q in the Plas- 
tic Ball spectrometer had an average value of a{q) ~ 10.5 
MeV/c and 12.5 MeV/c for the two-proton and two-deuteron 
systems, respectively. For small values of q, the resolution was 



around 7.5 MeV/c for the two-proton system. A bin size of 7.5 
MeV/c was chosen for the one-dimensional correlation anal- 
yses. It should be noted that the one-dimensional correlation 
functions presented here have been studied as a function of the 
momentum difference q rather than the more commonly used 
invariant four-momentum difference qi nv = \ \f—{Pi — P^) 2 - 
The experimental resolution for the pp (dd) system in qi nv at 
20 (40) MeV/c is about 15 (20) MeV/c. It increases to 20 (27) 
MeV/c at qi nv =100 (200) MeV/c and approximately levels off 
above that. 



3.1 Final state interaction 

To calculate the pair-wise final state interaction we used two 
computer codes for the proton-proton system: (i) the CRAB 
code written by S.Pratt ll27l which includes Coulomb and strong 
final state interactions, and (ii) the static model developed by 
A. Deloff 128I291 with both Coulomb and strong S-wave in- 
teractions included. The latter incorporates a potential of the 
delta-shell form 



2(j,V(r) =-(s/R)6(r-R), 



(2) 



characterized by the range R and the dimensionless parame- 
ter s representing the strength of the force (p is the reduced 
mass). Only the 1 Sq partial wave in the strong interaction was 
retained. This pp interaction is well known experimentally and 
the measured phase shifts can be satisfactorily reproduced up 
to q=150 MeV/c by taking the delta-shell potential parameters 
as s=0.906 and f£=1.84 fm. Both codes gave consistent results 
for the proton-proton system. 

For the deuteron-deuteron case the model by A. Deloff 1 28 ] 
was used in which the strong s-wave interaction in S=0,2 spin 
states was included together with the Coulomb repulsion be- 
tween two extended-size deuterons. The Coulomb interaction 
between two deuterons is different from that between two pro- 
tons because the deuteron is composite and relatively large. For 
two extended-size deuterons located at ri and r-2, the Coulomb 
potential can be written as 



V{\v 1 -v 2 \) =a / d 3 Xl d 3 x 2 



p(|xi-r 1 |)p(|x 2 -r 2 |) 
|xi— x 2 | 



(3) 



where p{x) is the charge density and a is the fine structure 
constant. Since the right hand side of Eq. (01 is translationally 
invariant, the potential V(r) depends only upon the difference 
r = |ri — r2 1, and for r larger than twice the deuteron radius 
takes the usual point-like form V(r) — a/r. For an assumed 
uniform charge distribution, this integral can be obtained in an- 
alytic form as: 



V(r) 



x[l- 



-3(l-x) 4 
i|(l-z)(5 



x)}} 0<a;<l (4) 
x > 1 



where x = r/Rc and Rc = 3.86 fm is twice the deuteron 
radius. A comparison of the Coulomb potential for extended- 
size and point-like deuterons is shown in Fig. [3] The Coulomb 
repulsion is observed to be substantially reduced when the ex- 
tended deuteron size is taken into account. 
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The one-dimensional two-particle correlation calculaE^fiTR')i 
based on a static spherical gaussian shaped source, specified by 
radius R : ( v j) 



p(x) = (27ri?2)-3/2 e ~x 2 /2Rl 



(5) 



In the case of a non-zero emission lifetime this approach gives (vii) 
an upper limit for the spatial extent of the source. 

To account for the Plastic Ball resolution all of the calcu- 
lated correlation functions presented in this paper were calcu- (viii) 
lated as: 



C(q)= / r(q,q')C(q')dq' 



(6) 



with the predicted correlation function C(q') convoluted with 
the Gaussian resolution function of the Plastic Ball r(q, q'). 



1 



] 1 (q-q') 

e 2 



') 2 



(2tt) 3 / 2 a(q) 
as determined from the GEANT simulations 

3.2 Systematic errors 



(7) 



Non-conservation of momentum in the background pair cor- 
relations: local (transverse) momentum conservation in real 
events, although not expected to be exact, may influence the 
correlation between real particle pairs due to the fact that a 
significant fraction, but not all, of the particles in the target 
rapidity region are measured in the Plastic Ball. To inves- 
tigate the sensitivity of the fitted parameters to the effect 
of momentum conservation, in addition to the requirement 
that the full background event have the same number of 
particles as the real event detected in the Plastic Ball, it was 
further required that the background event have a total mo- 
mentum of all measured particles equal to that of the real 
event: Pxot (measured) = Prot (reference) ± A, where A 
should be in the ideal case. In practice, A was chosen to 
be 300 MeV/c and 600 MeV/c for the two-proton and two- 
deuteron background events, respectively [0.07, 0.06 fm]. 
The granular structure of the Plastic Ball spectrometer: the 
GEANT program has been used to study the distortion of 
the correlation function due to the granularity of the Plastic 
Ball detector [0.15, 0.10 fm]. 

The bin width of the correlation function: to examine the 
dependence of the radius parameters on the q bin width, the 
analysis was also performed with bin widths of 10 MeV/c 
and 5 MeV/c [0.06, 0.06 fm]. 

Detector variation: the analysis was performed separately 
in two different azimuthal angle intervals: 0° — 180° and 
180° - 360° [0.01, 0.02 fm]. 

The effect of possible mis-identification of protons and deuterons 
due to the multiple hits to the same Plastic Pall module. 
Since the multi-hit probability depends very strongly on 6, 
the analysis was performed with two subsamples with 8 in- 
tervals of protons and deuterons: 60° - 90° and 90° - 160°. 
The ratio of contaminations due to multi-hits in these 6 in- 
tervals is about 10:1. The GEANT simulations also demon- 
strated that the contamination from spurious protons and 
deuterons contribute uniformly to their correlation func- 
tions. As a result, this background contamination does not 
affect the shape of the correlation function, leading to the 
same, within errors, radius parameters [0.15, 0.10 fm]. 



The total systematic uncertainty was calculated by adding 
in quadrature the contributions listed above. 



The systematic errors have been evaluated by comparison of 
the results obtained under different conditions. The systematic 
uncertainties are mainly introduced by the granularity of the 
detector and, to some extent, by the procedure used to construct 
the background pair distribution. The following sources of sys- 
tematic errors have been investigated leading to variations (rms 
estimates) of the radius parameters given in brackets, with the 
first and second numbers referring to the pp and dd system, 
respectively: 

(i) Different periods of data taking during the run [0.05, 0.05 
fm]. 

(ii) Different particle identification windows [0.09, 0.12 fm]. 

(iii) The width of the transverse energy Et bin used for cen- 
trality selection in the background pair distribution calcula- 
tions: the centrality bin size was halved and the result was 
compared to the nominal bin size result [0.01, 0.05 fm]. 



3.3 Two-proton system 

The measured two-proton correlation function is shown in Fig. [4] 
The solid line represents the result of the theoretical calcula- 
tion which gives the best fit to the data using the delta-shell 
potential. The final state Coulomb and strong interactions, the 
Fermi-Dirac effect, and the experimental resolution have been 
implemented in the theoretical calculation. The resulting radius 
parameter is 3.14± 0.03(stat.)± 0.21(syst.) fm. 

Although the delta-shell potential may not be realistic, its 
main advantage is that both the phase shift and the enhance- 
ment factor can be obtained in an analytic form. Furthermore, 
it is desirable to use the same model for calculation of both pro- 
ton and deuteron correlation functions. It was demonstrated in 
[29 1 that the calculated pp phase shifts as obtained from both 
the delta-shell potential and the REID soft core potential [ 30 1 



agree well with the experimental data, up to about 150 MeV/c. 
As a cross check, the REID potential was also used to fit the ex- 
perimental two-proton correlation function. The resulting cor- 
relation function is shown in Fig.@]as the dashed line. The fitted 
radius parameter is 3.05 ± 0.02(stat.)± 0.24(syst.) fm which, 
within errors, is consistent with the delta-shell potential calcu- 
lation. The calculated two-proton correlation functions without 
(SFulomb i nteraction and without any final state interactions are 
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Fig. 4. The measured two-proton correlation function. The errors 
showiQa£p statistical only. The, solid, and dashed lines represent the 
results of cQculatigtQ whi4£provi6£)the bBOfit tdllfltOlat^ kBsQig aJ40 
natively the delta-shell or REID potential, respectively. Thj 
correlation function for free propagation or with only 
interaction included are shown as dot-dashed and dotted lines. 



In order to make a direct comparison with other experimen- 
tal results, the CRAB package [27], provided by Pratt, was also 
used to extract the radius parameter for the two-proton sys- 
tem. The CRAB calculation takes into account the Coulomb 
and strong final state interactions. To describe the latter for 
the two-proton system the REID potential was taken. The mea- 
sured momentum distribution of protons was used as input to 
the CRAB calculation. To obtain a radius parameter, the cal- 
culation was performed for a set of different radius parameters 
and then compared to the experimental data by calculating the 
X 2 value. This gave a best fit radius parameter of 2.83 ± 0.20 
fm, that is consistent, within errors, with the delta-shell model 
calculation. The CRAB calculation result is shown in Fig. [5] 
The experimental correlation as a function of qi nv is also shown 
in Fig.|5]with the CRAB calculation. The solid and dashed lines 
represent the CRAB calculations with and without the experi- 
mental resolution effect included, respectively. 

The polar angle (or rapidity) dependence of the proton source 
radius parameter was also investigated. Within errors, the two- 
proton source radius, i?o showed no rapidity dependence, as 
seen from Fig. [6] over the entire polar angle interval. 

Two-proton correlations have also been studied in the tar- 
get fragmentation region at the SPS in the WA80 experiment 
using the Plastic Ball spectrometer. Measurements were made 
with 200 A GeV proton, Oxygen, and Sulphur projectiles, on 
targets of C, Al, Ag, and Au [11]. The source radius extracted 
from the present analysis of the two-proton correlation for cen- 
tral Pb+Pb collisions at 158 A GeV is compared to the radius 
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parameters extracted in the WA80 experiment in Fig. [7] The8 
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radius parameters are plotted as a function of A^ al t . All radii 
have been converted to root-mean-square radii. Within errors, 
the radius parameters measured in the WA80 and WA98 ex- 
periments at the SPS closely follow the target nuclear radius 

1 /3 

dependence of Aj, et with little apparent dependence on the 
projectile mass. Since the size of the participant overlap region 
depends on both target and projectile mass, the observation of 
little projectile mass dependence indicates that the protons de- 
tected in the target fragmentation region are emitted from the 
entire target volume, following substantial rescattering [11]. 



3.4 Two-deuteron system 

Figure [8] shows the measured two-deuteron correlation func- 
tion and the theoretical calculation which gives the best fit to 
the data. The solid and dashed line represent the correlation 
functions calculated with final state interaction taken into ac- 
count and that for free propagation. The fitted radius parame- 
ter is 1.59± 0.02(stat.)± 0.20(syst.) fm. The range R and the 
strength s for the delta-shell potential are found to be 0.52 ± 
0.01 fm and 1.08±0. 01, respectively. The measured two-deuteron 
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that since the correlation function C(q) — l + p(2q) 2 was 
for the calculations presented in section 3, a factor \[2 
should be applied to the radius parameters obtained from the 
one-dimensional analysis results presented above when com- 
parir g with the radius parameters obtained from the following 
mult -dimensional analysis. 

The multi-dimensional analysis has been performed in the 
Longitudinally CoMoving System (LCMS), in which the lon- 
gitucinal pair momentum vanishes (p 2 i + p Z 2 = 0). 

In order to extract that part of the correlation function of 
identical particles that is due to quantum statistical effects only, 
the measured multi-dimensional correlation functions have been 
corrected" for the effect of the strong and Coulomb final state 
interactions (FSI). The corrected correlation function is defined 

C(pi,V2) = K C raw (p 1 ,p 2 ) (9) 

whe7eK — Cfree(Q)/Cfsi(Q) is the correction factor, and 
Cfree(Q) and Cfsi{Q) are the calculated correlation func- 



protons have been measured in the target fragmentation regi^j with 
-0.6 < y <0.6 and 100 MeV/c < p T < 650 MeV/c . The sym- 731 ^ 6 * for free propagation and with final state interactions in 
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4.1 Two-proton system 



The Ilertsch-Pratt (BP) parametrization 



In 



tion 



the Bertsch-Pratt parameterization Q1I32I . the correla- 
fiinction for two fermions has the form: 



C = N ( 1 - \c-Qr R '-Q 2 °- R2 s-Qo- R °- 2 - 



(10) 



A ^rrett: N is the normalization constant and < A < 1 is the 
strength of the correlation which depends on the chaoticity of 
the emitter and on the purity of the measured data sample. The 
— strong fsi + coulomb for extended-size DEUTERgflS^ntum difference vector Q of two particles is decomposed 

into a longitudinal component along the beam axis (Qi) and a 
component transverse to the beam axis (Qt)- The transverse 
momentum difference is further decomposed into the outward 
and sideward components. The outward and sideward momen- 
tum-differences (Q ) and (Q s ) are defined as the components 



STRONG FSI ONLY 



FREE 



me|asured two-deuteron correlation function. The errors 
rrly. Th b .solid lin e 1 repres e nts 1 the ' r d sult .of 1 th e 



ides 



est 



the 



wit 



interaction included. The dashed and dot-dashed lines show 
relation function with only strong final state interaction taken 
count and that for free propagation, respectively. 



mem 



Q0final2flb 1c4Qie transverse momentum difference parallel and perpen- 
r to the total transverse momentum of the pair, respec- 
The sign convention is that Q s is always positive, Qi 
and Q are allowed to be positive or negative. R 2 t is the "out- 
longitudinal" cross term ll33l . Its value can be either positive 
or negative. In the LCMS system the average resolutions in the 
Q s , Qu and Q Q components for the pp system are 20, 23, and 
24 MeV/c, respectively. 

Figure|9]shows projections of the measured two-proton cor- 
relation functions, corrected for final state interactions (Eq. ©), 
onto the Qi, Q s , and Q Q axes. Also shown are the correspond- 
ing projections of the fit of Eq. (flOl . shown as solid lines. The 
projected results and the fit results for a given Q component 
have been integrated up to 75 MeV/c over the other two \Q\ 
mentum difference Q = (pi— P2) is used in the multi-dimensionabomponents. 

analysis. In this case the two-particle correlation C(Q) is given The parameters extracted from fits to the data using either 
by: the delta-shell corrected and CRAB corrected correlation func- 

C(Q) — 1 + p{Q) 2 ■ (8) tions are presented in Table [T] Also listed in this table are the 



correlation and extracted radius parameters show no depen- 
dence on the polar angle, similar to observations for the two- 
proton correlation (see Fig. 0. 



4 Multi-dimensional analysis 

For consistency of notation with previous analyses, the mo- 
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Table 1. Fitted parameters of the Bertsch-Pratt parametrization of the two-proton correlation function using the delta-shell model or CRAB 
model final state interaction calculations. Errors are statistical + systematic. 




a 
o 

Fig. 9QPgjections of the measured two-proton correlation function 
corrected wit h the (JRAB calculation fof the FSI. The solid lines show 



the projections of2£iGt with" tti06ertsch-pQtt paramelrQOion. Pi * 
jections for each Q component, for both the experimentffNdati 
the fit, have been integrated up to 75 MeV/c over the otnerOwo 
components. The errors shown are statistical only. 



'£ nf 4lk _ 

^ U Mine correspond! &R$&rA 

V/c) 



where (3t is the transverse velocity of the particle pair. Notice 
that Eq. (fTTt is valid only when opacity effects and transverse 
flow are ignored. The excess of R a over R s is due to the dura- 
tion of emission from the source in which particles with trans- 
verse velocity (3t travel a distance /3tAt towards the detector, 
resulting in an apparent extension of the source in that direc- 
tion. The sideward radius parameter R s is not affected by the 
duration of emission since it is perpendicular to 3t- 

Within errors, the radius parameters (Ri , R s , R a ) of the BP 
fit are quite close to each other (Table [TJ which suggests a vol- 
ume emission of the protons. The extracted R Q and R s param- 

...... ..... ,. _ .^eters, and the measured average value of 8t = 0.278 ± 0.09 

The statistical significance of the fit is demonstrated in Fig. 1101 r . . nr£ , r.u 

,. , , , j, t, „ — tor the proton pairs, gives 11.5 ± 2.5 fm/c as an estimate of the 

emission duration. 

In the case of a Bjorken scenario, with longitudinal expan- 
sion of the system without transverse flow, the wit dependence 
of the longitudinal radius reflects the lifetime of the source[36|: 



fitted parameters without the cross term R 2 ol . Within errors, the 
delta-shell and CRAB corrections lead to the same values of fit- 
ted parameters. The fitted radius parameters with and without 
the cross term do not differ within errors. The cross term R a i 
vanishes in the LCMS frame for longitudinally boost-invariant 
systems. The observation that the cross-term R? ol is small, in- 
dicates that the proton source undergoes an approximate boost 
invariant expansion. 



where contour lines of equal confidence level for R a versus R 
are plotted. The spacing of contours corresponds to one stan- 
dard deviation in x 2 value. The solid line, which represents 
equal sideward and outward radii, is 9 standard deviations from 
the minimum \ 2 value. 

For transparent sources, where particles are emitted from 
throughout the entire volume, the longitudinal radius Ri is in- 
terpreted to measure the longitudinal size of the system di- 
rectly, while the sideward radius provides a measurement of the 
transverse size. The duration of emission At is related B34I35I 
to the sideward radius R s and outward radius R a as: 



R l = T f( T f/ m T) 



1/2 



(12) 



1 



Ar = —^/Rl-Rl 

PT 



(ID 



where Tf is the freeze-out temperature and ttit — \l m? + p T . 
Under the assumption that Tj = 120 MeV, a source lifetime of 
about 13 fm/c is obtained. 

Strictly speaking, the interpretation of r/ as the freeze-out 
time is based on the assumption that the longitudinal boost- 
invariant velocity profile existed not only at the time of freeze- 
out but also throughout the dynamical evolution of the reaction 
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zone. If this were not the case, then r/ would be larger. Hence, 
Tf estimated by Eq. (fT2l > gives a lower limit for the lifetime of 
the source HI. 

The Yano-Koonin-Podgoretsky (YKP) parametrization 



In the Yano-Koonin-Podgoretskvparameterizafion ll37l38l . 
the two-proton correlation function is defined by the formula: 

C = N (l — \e-^T R T~(Qi-Q 2 o)^-(Q- u ) 2 ( R2 a+^^ (13) 

where U = 7(1, 0, 0, (3) is a 4-velocity with only a longitudinal 

1/ 



component, 7 = — {3 2 , Qt and Qi are the components 

of the momentum difference projected onto the transverse and 
longitudinal directions, respectively, Qo is the difference in en- 
ergies. The sign convention is Qo is always positive and Qi can 
be either positive or negative. This parametrization has an ad- 
vantage in that the three YKP radius parameters do not depend 
on the longitudinal velocity of the measurement frame, and the 
velocity (3 is closely related to the velocity of the effective par- 
ticle emitter. 

The parameters extracted from YKP fits to the delta-shell 
corrected or CRAB corrected correlation functions are listed in 
Table |2] Within errors, the fitted parameters obtained with the 
two different final state interaction corrections are consistent 
with each other. 



Table 2. Fitted parameters of the YKP parametrization of the two- 
proton correlation functions using either the delta-shell model or 
CRAB model FSI corrections. Errors are statistical + systematic. 





Delta-shell Model 


CRAB Model 


A 


0.75±0.04±0.04 


0.77±0.03±0.03 


Ri(fm) 


4.64±0.16±0.36 


4.97±0.08±0.46 


R T (fm) 


4.42±0.13±0.40 


4.73±0.06±0.46 


fti(fm) 


10.54±0.46±1.18 


10.58±0.46±1.12 


13 


0.032±0.023±0.010 


0.036±0.014±0.011 


X7ndf 


1.42 


1.47 



4.2 Two-deuteron system 



The Gaussian parametrization in the LCMS 

Fig.QT]shows the experimental one-dimensional two-deuteron 
correlation function in the LCMS frame, after correction for fi- 
nal state interactions using the delta-shell model described in 
section 3. The correlation is fitted to the Gaussian correlation 
function: 



C(Q) = N(1 + \c-Q 2r1 



(14) 



The fitted parameters are: A = 0.43±0.03(stat.)±0.08(syst.) 
R = 2.50 ± 0.10(stat.)±0.28(syst.) fm. The radius parameter is 
(njagreerr ent within errors with the radius (R = 1.59 x \[2 
£25. fm) extracted by fit to the M§±aMfeM-&tfkf^SHfelatio|i 
function with the full delta-sheH mad@i(itall<atbtt}i9i2(8^f^g.|8) 




0.8 



Fig. QlPrhe fitted two-deuteron correlation function as a function of 
momentun l difference Q in the LCMS frame after correction for fina 1 
state intergptions. 50 1 QO 150 200 250 300 



Q(MeV/c) 



The Bertsch-Pratt parametrization 



Since the YKP velocity j3 is close to 0, it demonstrates that 
the LCMS frame coincides with the YK frame (the frame for 
which the YKP velocity parameter j3 vanishes). For transpar- 
ent sources, and in the absence of flow, the radius parameters 
in the YK frame have a convenient physical interpretation. Ri 
and Rt are then interpreted as measures of the longitudinal 
and transverse size of the source, while Rq reflects directly the 
duration of emission from the source. 

With similar fitted Ri and Rt parameters, and longitudinal 
velocity [3 compatible with 0, the YKP fit indicates a symmetric 
emission geometry with approximate boost invariant expansion 
of the source. The emission duration extracted from the YKP 
or BP fits are in agreement within errors. 

The fitted YKP parameters are seen to be consistent with 
the BP parameters. Since the BP and YKP parametrizations 
are mathematically equivalent and differ only in the choice of 
independent components of q, the agreement between the two 
sets of fitted parameters serves additionally as a consistency 
check of the fitting procedures. 



The two-deuteron correlation function has been corrected 
by the calculated delta-shell model two-deuteron final state in- 
teraction and fitted with the two-boson correlation function us- 
ing the Bertsch-Pratt parametrization: 

C = N (l + \ e ~Qr R2 -Q 2 - R2 -Q 2 -Q 2 - 2 -Q°Q'- R2 i^ . (15) 

Projections of the measured two-deuteron correlation func- 
tion on the Qi, Q s , and Q axes are shown in Fig. [T2] The 
measured and fitted projections in a given Q component have 
been integrated up to 75 MeV/c in the other two \Q\ compo- 
nents. All Q components have approximately the same average 
resolution of 32 MeV/c. 

The fitted parameters, with and without the cross term R?,, 
are listed in Table [3] Within errors, the values of parameters 
extracted from the two different fits are the same. Similar to 
the result for the proton source, the small cross-term indicates 
that the deuteron source also experienced an approximate boost 
invariant expansion. 
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A 


0.29±0.02±0.09 


Ri(fm) 


1.29±0.07±0.42 


R T (fm) 


2.35±0.05±0.39 




-46.95±1.69±8.82 


P 


0.020±0.002±0.007 


X Vndf 


1.42 
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Fig. 12. Projections of the measured two-deuteron correlation func- 
tion dOrSe ted 1 with the 1 dclta- ' sh c ll rhodcl ' calculation for thqlFSI ahd 



Q, l< 75 MeV/c 

» , _,, . .. ' again indicates that the deuteron source undergoes an approxi- 
Q, < 75 MeV/c? . . 6 FF 

' mat£ boost invariant expansion. 



the YK frame the YKP radius parameter Rq can be ex- 



press ed as (41 j: 



fitted wittiT&PQertsch-Pratt' pJ?ametrizationQ Projections ftPQach Q 200 
component for both the experimental data and the fit hav^eenf f^ff^^ /q\ 
grated over the range 0-75 MeV/c in the other two \Q\ compOftkents. ' 
The errors shown are statistical only. 



PI 



(R - R s 



(17) 



Table 3. Fitted parameters of the Bertsch-Pratt parametrization of the 
two-deuteron correlation functions using the delta-shell model FSI 
correction. Errors are statistical + systematic. 





with cross term 


without cross term 


A 


0.37±0.02±0.19 


0.37±0.02±0.19 


Ri(fm) 


2.17±0.10±0.24 


2.13±0.08±0.18 


R,(fm) 


3.14±0.11±0.78 


3.10±0.10±0.67 


Roifm) 


0.01±0.08±0.02 


0.01±0.04±0.03 


Roi(fn?) 


0.20± 0.08 ±0.69 




XVndf 


1.52 


1.52 



Within errors, the fitted Ri and R s parameters are close to 
the radius value obtained from the one-dimensional fit, while 
R is consistent with 0. For a transparent, azimuthally sym- 
metric emission source, it is expected that R s and R Q should 
be similar, except that R will be extended by the duration of 
emission, as discussed with respect to Eq. ( fTTT i. The observa- 
tion that R << R s is inconsistent with such a transparent 
source, but may result naturally from surface emission from 
an opaque source |39 40]. Since deuterons are relatively large 
weakly bound objects, it might be expected that only those 
deuterons produced on the freeze-out surface via coalescence 
of an emitted proton and neutron survive. This would naturally 
lead to an opaque emission source, as observed. 

The Yano-Koonin-Podgoretsky parametrization 

The two-boson correlation function in the Yano-Koonin- 
Podgoretsky parametrization is given by: 

C = N (l + \ c -Q 2 tK 2 t-(Qi~Q 2 )Ri-(Q-u) 2 (R 2 b +R^ _ (16 ) 



where /3± is the velocity of the particle pair transverse to the 
beam direction, and R Q and R s are the outward and sideward 
radii of the BP parametrization. The negative value of Rq w 
—47 fm 2 obtained from the YKP fit is in agreement with the 
observation of strong opacity 1 1 39 42 1 of the deuteron source 
made in the previous section. 



5 Discussion 

Although one does not necessarily expect boost-invariance in 
the target fragmetnation region at SPS energies, the small val- 
ues of the cross-term R i in the BP parametrization as well as 
the consistency with zero longitudinal velocity (3 in the YKP 
parametrization of the correlation functions, indicates that both 
proton and deuteron sources exhibit boost-invariant expansion. 

The measured one-dimensional radius parameter extracted 
for protons is markedly larger than that for deuterons. Compar- 
ison of the present results for protons and deuterons with those 
for 7r + , also measured in the Plastic Ball |43 1, and 7r~ measured 
in the WA98 negative-charged particle spectrometer [4445], 
and with K + measured in the NA44 experiment [46 1 for cen- 
tral Pb+Pb collisions at 158 A GeV, reveals that the radii are 
observed so show a mass ordering with R^^ > Rkk > Rpp > 
i?dd- This observation is in agreement with previous reports 
[47 1 that lighter particles tend to give larger radius parameters. 
It is interesting that a large loosely bound composite object like 
the deuterons follows the same trend as observed for elemen- 
tary particles. It has been shown by Alexander [4-7 48 1 that the 
source-size mass dependence of hadrons emerging from Z° de- 
cays produced in e + e~ annihilation can be reproduced using 
the Heisenberg uncertainty relations to derive the relations: 



R 



A 



rn 



(18) 



The parameters of the YKP fit to the two-deuteron correla- 
tion are listed in Table|4] The longitudinal velocity (3 is close to 
indicating that the LCMS coincides with the YK frame and 



and 



(19) 
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where A — c^hAt is a time scale constant, Ri is the longitudi- 
nal radius parameter, and m and tot are the mass and the mean 
transverse mass, respectively. An alternative explanation using 
a QCD derived potential B71 proved to be equally successful. 
A different approach to explain the R and Ri mass dependence 
was given by Bialas and Zalewski [49 50 5 1 1. In this descrip- 
tion the radius parameter of the source is mass independent 
and its apparent decrease is a consequence of the momentum- 
position correlation expressed in the Bjorken-Gottfried condi- 
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led to the conclusion that this 
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Fig. l£ The rriT dependence of the longitudinal radius parameter in 



central ^P b+Pb collisipns at 158 ^ GeV. Thq error bars indicate the 
statisticalQand systematic errors adlded in quadrature. The sslid curve B.5 




:ta dc 



the result of a fit of R — A/ 

with the Plastic Ball spectrometer, while the dashdd'3u*Ve n ThcfiMe: 
also mid-rapidity n~ data from WA98 and kaon data from the NA44 
and NA49 experiments. 




Fig. 13: The dependence of radius parameters extracted 



dimensio nal two-particle correlations on the particle ma^s in central 
Pb+Pb tojlisions at \p§ A GeV TJie error barsirepresent tjie statistical 5 
and systematic errors added in quadrature. The solid curve, is^the resujt 
of fit of Eq. ([Till to the three WA98 data points measiflidlwiffi PKfitid 
Ball spectrometer. The dashed curve includes also the n~ data from 
WA98 and kaon data from the NA44 experiment. 



data are compared to fits of the form Ri = Aij ^JrriT, where 
i = l,T. As can be seen from Figs. [14] and Q3] within exper- 
imental uncertainties, the fitted Ri and Rt radius parameters 
are consistent with A/^/ttit scaling. A fit to the three WA98 
data points measured near mid-rapidity with the Plastic Ball 
sgectrdn ^Her yields A t = 3.05 ± 0.17 fin GeV 1/2 with x 2 /ndf 

^= 1 8 



rr 




Radius parameters extracted from one-dimensional two-particle 
correlation functions for central Pb+Pb collisions at 158 A GeV 
are shown as a function of the particle mass in Fig. Qj] The 
7r~ 7r~ radius parameter was measured with the negative-charged 
particle spectrometer of the WA98 experiment [44 45 ] and the 
two-kaon radius parameter was reported by the NA44 Collab- 
oration 11461 . both measured near mid-rapidity (y) ~ 2.9. In 
general, the radii follow the y/m dependence expectation of 
Eq. ([Ull. A fit of Eq. (Q~8]i to the 7r+ ED, proton, and deuteron 
radii measured with the WA98 Plastic Ball gives A — 2.84 ± 
0.13 fm GeV 1 / 2 with x 2 /ndf = 12, as shown by the solid curve. 
This corresponds to a time scale of At — 40.8 fm/c using the 
above expression for A, or to a freeze-out lifetime of r/ = 8.2 
fm/c using Eq. ([1211 . A fit to all data points (dashed line) results 
in A = 2.75 ± 0.04 fm GeV 1 / 2 with x 2 /ndf = 7.2. 

Previous results on pion and kaon interferometry for S+Pb 
and Pb+Pb reactions at the SPS have demonstrated scaling of 
the Ri and Rt radius parameters on the transverse mass, tot 155156 
The longitudinal and transverse radius parameters from the multi- 
dimensional YKP analysis of the two-proton and two-deuteron 
correlation functions are plotted as a function of the mean trans- 
verse mass in Figs.[l4]and[T5] Also shown are radius parame- 
ters for 7r + measured with the Plastic Ball l43l . Radii measured 
near mid-rapidity for ir~ in WA98 M44I45L and for kaons from 
the NA44 [56| and NA49 [57] experiments are also shown. The 



all data points are incluied;-^! %5^g the fit yields 
0.05 fm GeV 1 / 2 with ^ 2 /ft(^^M98 
pp, WA98 



dd, WA98 
K + K + , NA44 
K + K + , NA49 
K"K", NA49 




The m,T dependence of the transverse radius parameter- -in 
Pb+Pb collisipns at 158 ^ GeV. Thq error bars pdicate the 



Fig. ii 

central 

statistic'alOand systematic errors atlded in quadrature. The2solid curvg.5 
is a fit of Rt = A/ \frnr to the WA98 data points measured with tl 
Plastic Ball spectrometer only, while the dashed 
mid-rapidity 7r~ data from WA98 and kaon data from the NA44 and 
NA49 experiments. 



Although there is no theoretical justification to expect an 
Rt=A/ \jm,T dependence, the data are seen (Fig. [15]) to follow 
such a dependence, within errors. A fit to the three WA98 data 
points measured near mid-rapidity with the Plastic Ball spec- 
trometer yields A T = 3.38 ± 0.18 fm GeV 1 / 2 with x 2 /ndf = 
2.4 When all data points are included in the fit, the fit yields 
A T = 2.68 ± 0.04 fm GeV 1 / 2 with x 2 /ndf = 3.3. A fit to the 
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pion and kaon data from NA44 experiment alone gave At = 
3.0 ± 0.2 fm GeV 1/2 (56). 



6 Summary 

In summary, we have measured two-proton and two-deuteron 
correlation functions in the target fragmentation region for cen- 
tral Pb+Pb collisions at 158 A GeV. 

In the one-dimensional analysis, the radius parameters were 
extracted from the measured correlations using calculations that 
assume a static Gaussian-shaped source. The extracted proton 
source radius parameter is about two times larger than that of 
the deuteron source. Comparison of the proton and deuteron 
radius parameters with those extracted from one-dimensional 
two-pion and two-kaon correlations reveals a common A/y/m 
dependence with A « 3 fm GeV 1 / 2 . Such a dependence can be 
explained as a consequence of the Heisenberg uncertainty rela- 
tions, under the assumption of a common duration of emission. 

The multi-dimensional analysis demonstrates that the pro- 
ton source exhibits a volume emission with long emission time, 
whereas the deuteron source is strongly opaque. The cross-term 
R 2 ol from the Bertsch-Pratt parametrization fit and the longitu- 
dinal velocity (3 from the Yano-Koonin-Podgoretsky parametriza- 
tion fit are consistent with zero, suggesting that the proton and 
deuteron sources undergo an approximate longitudinal boost 
invariant expansion. 

The longitudinal and transverse radius parameters, Ri and 
Rt, extracted from the multi-dimensional correlation analysis 
follow a common Af^/mr scaling for pions, kaons, protons, 
and deuterons, with A w 3 fm GeV 1 / 2 in both cases. The exis- 
tence of a universal function describing the dependence of the 
radius parameters on the transverse mass for different particle 
species in the mass interval from pion to deuteron may indi- 
cate an approximately simultaneous freeze-out of the studied 
hadrons. 
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